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Abstract 



A model of localized classical electrons coupled to lattice degrees of free- 
dom and, via the Coulomb interaction, to each other, has been studied to 
gain insight into the charge and orbital ordering observed in lightly doped 
manganese perovskites. Expressions are obtained for the minimum energy 
and ionic displacements caused by given hole and electron orbital configura- 
tions. The expressions are analyzed for several hole configurations, including 
that experimentally observed by Yamada et al. in Lay/gSrx/gMnOs. We find 
that, although the preferred charge and orbital ordering depend sensitively on 
parameters, there are ranges of the parameters in which the experimentally 
observed hole configuration has the lowest energy. For these parameter values 
we also find that the energy differences between different hole configurations 
are on the order of the observed charge ordering transition temperature. The 
effects of additional strains are also studied. Some results for La^Ca^MnOa 
are presented, although our model may not adequately describe this material 
because the high temperature phase is metallic. 

71.38.+i, 71.45.Lr, 71.20.Be, 72.15.Gd 
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Over the last few years much attention has been focused on manganese perovskite-based 
oxides, most notably the pseudocubic materials Rei_ x Ak x Mn0 3 . (Here Re is a rare earth 
element such as La, and Ah is a divalent alkali metal element such as Ca or Sr.) The 
initial motivation came from the observation that for some range of x, and temperature, T, 
resistance can be reduced by a factor of up to 10 7 in the presence of a magnetic field.0 Two 
other interesting physical phenomena occurring in this class of materials are charge ordering 
and orbital ordering.! In this paper, we study the connection between the two. 

The important electrons in Re\^ x Ak x yhi0 3 are the Mn e g electrons; their concentration 
is 1 — x. For many choices of Re, Ak, and x, especially at commensurate x values, the e g 
charge distribution is not uniform and it indeed appears that a fraction x of Mn ions have 
no e g electron while 1 — x have a localized e g electron. A periodic pattern of filled and empty 
sites is said to exhibit charge ordering. There are two e g orbitals per Mn ion. A localized Mn 
e g electron will be in one linear combination of these; a periodic pattern of orbital occupancy 
is said to exhibit orbital ordering. Recently, Murakami et a/.i observed the charge ordering 
transition accompanying simultaneous orbital ordering in Lai/2Sr 3 / 2 Mn04 at 217 K (well 
above the magnetic phase transition temperature 110 K). It indicates that the interplay of 
the charge and orbital ordering to minimize the lattice energy could be the origin of the 
charge ordering. In this paper we present an expression for the coupling between charge and 
orbital ordering, with different charge ordering patterns favoring different orbital orderings. 
We also argue that the orbital ordering energy differences determine the observed charge 
ordering in lightly doped manganites. Localized charges induce local lattice distortions, 
which must be accommodated into the global crystal structure; the energy cost of this 
accommodation is different for different charge ordering patterns. 

To model the charge and orbital ordering, we assume that the electrons are localized 
classical objects, so that each Mn site is occupied by zero or one e g electron, and each e g 
electron is in a definite orbital state. This assumption seems reasonable in the lightly doped 
materials such as La7/sSri/ 8 Mn03, which are strongly insulating at all temperatures,!! but 
may not be reasonable for the Lai^Ca^MnOs composition]! where the charge ordered state 
emerges at a low temperature from a metallic state. We proceed by calculating the energies 
of different charge ordering patterns, emphasizing the 1/8 doping case. It is practically im- 
possible to consider all possible charge ordering configurations. Therefore, we consider the 
three configurations shown in Fig. which are the only ones consistent with the following 
basic features of the hole-lattice implied by the experimental results by Yamada et al% : 
invariance under translation by two lattice constants in the x or y direction, four in the z 
direction, and an alternating pattern of occupied and empty planes along z direction. The 
configuration in Fig. |I|(b) is the one proposed by Yamada et a/.i to explain their experimen- 
tal results for La7/gSri/gMn03. For localized electrons there are three energy terms : the 
coupling to the lattice, which will be discussed at length below, the Coulomb interaction, 
and the magnetic interaction. 

First, we argue that the Coulomb energy cannot explain the observed ordering pattern 
or transition temperature. We take as reference the state with one e g electron per Mn and 
denote by Sqi the charge of a hole on a Mn site. From the classical Coulomb energy 

1 x . 5qMj m 

(-^Coulomb — ~ — 2^ ' V 1 / 

Ze v£j r iJ 
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one finds that the difference in energy between the configurations in Fig. [I] is 



^<Coulomb, per hole — ~ — ' V / 

where o is a site containing a hole and A(<5g) is the difference in charge between the two 
configurations. We estimated the above infinite sum by repeated numerical calculations for 
larger and larger volumes of the unit cells around the origin. We find that Fig. 0(c) has the 
lowest energy; 12 meV/e D lower than Fig. |l|(b), and 27 meV/e D lower than Fig. |](a). 

To estimate the magnitude of the Coulomb energy differences, we need an estimate for the 
dielectric constant eo, which we obtain from the measured reflectivity for Lao.gSro.iMnOs,! 
and the Lyddane-Sachs- Teller relation! uo 2 L = u^eo/too- At frequencies greater than the 
greatest phonon frequency the reflectivity is close to 0.1, implying ~ 3.4; the reflectivity 
is near unity between ujt = 0.020 eV, and ujl = 0.024 eV, implying e ~ 5.0. Because both 
Lay/gSrx/gMnOs and Lao.gSro.iMnOs are insulating and have similar compositions, their 
static dielectric constants are expected to be similar. Using e ~ 5.0, the energy difference 
between different configurations of holes is only around 2.4 meV, or 30 K per hole, which 
is small compared to the observed charge ordering temperature of 150 K — 200 K of these 
materials. The inconsistency with the experimentally observed hole configuration and the 
smallness of the energy difference scale indicate that the electrostatic energy is not the main 
origin of charge ordering for this material. 

Even though the magnetic and charge ordering transitions show a correlation in 
Lay/gSri/gMnOa,! we do not think that the magnetic contribution to charge and orbital 
ordering is as important as the lattice contribution for three reasons. First, in undoped 
LaMnOs, the orbital ordering and the structural phase transition occur at around 800 K 
and the magnetic ordering at around 140 K,&0 suggesting that the magnetic effects are 
relatively weak. Second, in Lay/gSr^gMnOa the Mn spins are ferromagnetically ordered 
with moment close to the full Mn moment at temperatures greater than the charge or- 
dering temperature,! and ferromagnetic order does not favor one charge configuration over 
another. Third, although in La7/gSri/gMn03 antiferromagnetic order appears at the charge 
ordering transition, the antiferromagnetic moment is very small (less than 0.1 of the full Mn 
moment),! so the energy associated with this ordering must be much less than 140 K/site 
associated with magnetic ordering in LaMnOa. Therefore, we think that the canted anti- 
ferromagnetism occuring upon charge ordering in Lay/gSrx/gMnOa (Ref. 7) is not the cause 
but the effect of the charge and orbital ordering. We now turn our attention to the lattice 
energy. 

A classical model for the lattice distortions of the insulating perovskite manganites has 
been derived in Ref. 10, and shown to be consistent with experimental results on LaMnOa. 
This model is adopted here with an additional term, an energy cost for shear strain. We now 
briefly outline the model, which is explained in more detail in Ref. 10 and the Appendix. 
The ionic displacements included are the vector displacement 5i of the Mn ion on site i, 
and the a directional scalar displacement (a = x, y, and z) of the O ion which sits 
between the Mn ion on site i and the Mn ion on site i + a. For convenience, 5i and it" 
are defined to be dimensionless in the following way: the lattice constant of the ideal cubic 
perovskite is b, the Mn ion position in the ideal cubic perovskite is Ri, the actual Mn ion 
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position is Ri + b5i, and the actual O ion position is Ri + (b/2 + bu^)d. The lattice energy- 
is taken to be harmonic and depends only on the nearest neighbor Mn-0 distance and the 
first and second nearest neighbor Mn-Mn distances. The spring constants corresponding to 
these displacements are K±, K 2 , and K 3 as shown in Fig. [2|. Because K\ and K 2 involve 
bond stretching, while K 3 involves bond bending, K\ > K 2 ^> K 3 is expected. Thus, 

-^lattice = E Mn _ Q + ^Mn-Mn,first + -E'Mn-Mn.second > where 



E 



Mn-O 



E 



Mn-Mn, first 



E, 



Mn-Mn, second 




(3) 
(4) 



(5) 



In the above equations a denotes x, y, and z, and (a, b) represents (x,y), (y,z), and (z, x). 
-^Mn-Mn,second was not considered in Ref. 10. The shear modulus produced by this term is 
important, because without it, a Mn ion on site i + x can have arbitrary large y directional 
displacement relative to the Mn ion on site % at no cost of energy. For this reason, the 
model with K 3 = has singularities, whose proper treatment requires K 3 7^ in our model. 
However, still we expect K 3 will be much smaller than K\ or K 2 . Therefore, in order 
to simplify the calculation, the K 3 /Kx — > limit has been taken after the expression of 
minimized energy and equilibrium ionic displacements have been obtained. 

Second, we consider the electronic degree of freedom. We parameterize the electron 
density by the variable h{. If an electron is present on site i, hi = 0; if no electron is 
present, hi — 1. If there is an electron on site i, the electron orbital state, which is a linear 
combination of the two e g orbitals, is parameterized by an angle 9i as follows. 



|?M#j) >= cos6 , j|c?3 2 2_ r .2 > + sin 6i\d x 2_ y 2 > 



(6) 



with < $i < 71. The electron orbital state couples to the distortion of the surrounding 
oxygen octahedra through the Jahn- Teller distortion. The coupling is given by 



£jt = -A£(l - ^)[cos2^{< - ~« + vfi} + sin2^« 
= -AX)(l-/i i Kcos2(0 i + Va), 

i,a 



«?)] 



where 



-7r/3, IPy 



7T/3, ip 



0. 



(7) 



(8) 
(9) 



If a hole is present on site i, it attracts the surrounding oxygens equally, giving rise to a 
breathing distortion energy given by 



4 



E holc = /3A E fn{vf + V? + <). (10) 

i 

The parameter (3 represents the strength of the breathing distortion relative to the Jahn- 
Teller distortion. Finally, following Kanamori,! we include a phenomenological cubic anhar- 
monicity term given by 

^anharm = ~^E(1 - h { ) COS 6^. (11) 

i 

The sign has been chosen so that the electron orbital states of |3x 2 — r 2 >, |3y 2 — r 2 >, or 
|3z 2 — r 2 >, with x, y, and z pointing toward nearest oxygen ions are favored when A is 
positive. The total energy, which is the sum of all the above energy terms, is given by 

Etot — Eyin-O + -E-Mn-Mn,first + -^Mn -Mn, second + Ejt + Eh i e + -Eanharm- (12) 

We minimized E to t about c\ a 's and m-'s for fixed hole and orbital configurations. These 
are conveniently expressed in terms of the variables 52, u% hj;, and c| defined in the following 
way. 

<^ = X>-*%2, (13) 

k 

< = E e_i ^ u I» ( 14 ) 

k 

hi^e-t-Hg, (15) 

k 

(l-/l i )cO8 2(0 i + Va) =E e "^4 ( 16 ) 

k 

The details are shown in the Appendix. The minimized energy per Mn ion may be written 

as 



Etot £7 v r< a -^anharm 



where 



A 2 



+ K 2 (l - cos* a )K/% - <$)(ph_ n ~ C\), if k a ± 



fc 10, iffc a = 1 ; 



l2 



= -^T2^[ 3 (^o) 2 + E(^) 2 ]- (19) 

The long wave length strain e ab , and the 0) components of the ionic displacements are 
given as 

2A 

* b = - Kl + 2 K 3 {J3h ° (20) 

( A[Xi+K 2 (l-cosfc a )] l- e -'fca mi a \ :f J, / n 

7 ,2 —J (ifi+2i<r 2 )Xl 1-cosfea^ *i a ^ foil 

Vo " \ o, if fc a = 1 ] 

S& - I ~ % (A' 1+ 2X 2 Kl-cos k~) (P h k ~ C g)' lf ^^° (oo\ 

"lo, if k a = o. 
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Because h^s and (1 — hi) cos 2(6*; + ^a)'s are bounded by ±1, we cannot treat h%s and 
c|'s as independent variables to minimize E tot . Therefore, we minimize E tot over the orbital 
variable 0j's at fixed hole configurations; the ground state is then the hole configuration of 
the lowest energy. 

For La7/gSri/gMn03, we consider the three hole configurations shown in Fig. [I], each of 
which is a Bravais lattice, with a unit cell containing one Mn site with a hole and seven 
Mn sites without holes. The orbital configuration may be different in different unit cells 
of the lattice defined by the holes. We consider the case where the orbital configuration is 
the same in each unit cell. In addition to that, we also consider all possible two sublattice 
symmetry breakings. Therefore, we have seven (if no symmetry breaking) or fourteen (if 
two-sublattice symmetry breaking) orbital variable 9i's. E tot /N in Eq. ( |17|) for each con- 
figuration is expressed in terms of those variables through Eqs. (|TTD, (|15D, (|16|), (jl8D, and 
(|l9l), and is minimized about 0j's. For this minimization, we use the FindMinimum routine 
in Mathematica in the following way: for each set of parameters, and for each configuration, 
we check the local minimal values by using 50 — 200 random starting values of O^s. 

According to Ref. 10, X/Ki ranges over 0.04 — 0.05, and KijK\ is between and 1. 
A/Ki ranges around 0.0002, and K\ « 200 eV.0 Recently, a local breathing distortion 
of 0.12 A has been directly observed in Lao.75Cao.25MnO3.EJJ The Jahn- Teller distortion is 
estimated around 0.15 A from the Mn-0 distances of LaMn03.@ This implies that the 
breathing distortion and the Jahn- Teller distortion in these materials have similar order of 
magnitude, i.e., (3 = 0(1). We varied f3 in the range of — 10, and A/K% in the range of 
- 0.00035, with \jK x = 0.045, K 2 / K x = 0.5, and K x = 200 eV. For each set of those 
parameters, the minimum energy per hole for each fixed hole configuration in Fig. [1| has 
been found. By comparing them, we find the most favored hole configuration for each (3 and 
A/Kx, which is shown in Fig. [3] as a plot in (3-A/Ki plane. 

At large (3 (> 7), the configuration shown in Fig. |I](c) is the most favored, and that 
shown in Fig. [TJ(a) is the least favored. This can be related to the fact that, in y-z and z-x 
directional planes, Fig. |l|(c) has the most even distribution of holes, and Fig. 0(a) has the 
least even distribution. For large (3, the contraction of oxygen octahedra toward holes is 
strong, and an uneven distribution of holes generates larger strains and elevates minimum 
energies. Particularly, the square hole net squeezes the electron orbital at the center of the 
square along the direction perpendicular to the square plane. In the cubic hole configuration 
of Fig. |I|(a), the six squeezed electron orbitals point toward the cubic center, putting the 
electron orbital at the center at high energy, which is consistent with our result that Fig. 
|I|(a) has far higher minimum energies than Figs. |](b) and |](c) in large (3 limit. 

As f3 is decreased into the range of 2 — 5, the favored hole configuration becomes that 
of Fig. |l|(b), which is the experimentally observed hole configuration. We expect that the 
difference of the energy per hole between the ground state hole configuration and the next 
lowest energy hole configuration corresponds approximately to the charge ordering temper- 
ature. The calculation results indicate that, when f3 is in the range of 2.0 — 2.5 or around 
5.0 and AjK\ = 0.0002, the charge ordering temperature is around 100 — 200 K, which 
is consistent with experimental results. As (3 is decreased further, the most favored hole 
configuration changes further and the temperature difference scale decreases. 

Figure |3| also shows the tendency that the configuration of Fig. |](c) becomes more favored 
as A/Ki increases. We think this occurs because the anharmonicity energy distorts the 
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oxygen octahedra tetragonally, which can be more easily accommodated by the tetragonal 
hole configuration of Fig. 0(c). 

In Table Q, we have shown an example of the orbital states, ionic displacements, and 
uniform strains corresponding to the minimum energy configuration for Fig. |l|(b) when 
AjK x = 0.0002, X/Kx = 0.045, K 2 jK x = 0.5 and (3 = 2.5. x, y, and z directions are shown 
in Fig. |I|. The nearest Mn-Mn distance is unit, (nf , n\ , nf) is defined in such a way that 
(nf , nl nf) + ^(2, 0, 0) + 2iV 2 (0, 2, 0) + 2iV 3 (l, 0, 2)'s and (nf , nf , nf) + N x {2, 0, 0) + (2N 2 + 
1)(0, 2, 0) + (2N 3 + 1)(1, 0, 2)'s, where N x , N 2 , and N 3 are integers, represent the coordinates 
of the sites indexed by i. k = parts of the ionic displacements have been subtracted to 
find the non-uniform parts of the displacements. 

The energy expressions in Eqs. flI8|) and ( P"§|) are adequate for bulk materials. When 
the material is grown on a substrate as a thin film, generally there is a strain generated 
by lattice mismatch between the film and the substrate materials. To see the effect of this 
strain, we add a term proportional to Cq (a' = x, y, or z) to the energy, which corresponds 
to an a' directional strain. Using a parameter g, we replace £^, =Q in Eq. ( |i~9"D by the following 
expression: 

^=o = - Kl + 2K 2 Wk ° )2 + 9C ° + ? (C ° )2] - (23) 

We repeated similar calculations to find the favored hole configurations for different values 
of the applied strain, parameterized by g. The applied strain breaks cubic symmetry. Some 
of the hole configurations also break cubic symmetry. For these cases the energy depends on 
the relative orientation of the strain and hole symmetry breakings. We consider all possible 
orientations and find the lowest energy state. We have varied g between —0.4 and 0.4, and 
P between and 7, with AjK x = 0.0002, \jK x = 0.045, and K 2 /K x = 0.5. The results 
are shown as a phase diagram in (3-g plane in Fig. |4]. It shows that Fig. [TJ(c) configuration 
is favored more as \g\ increases. This feature can be understood in the following way. 
For small g's, the leading correction to the minimum energy for each hole configuration is 
— \ 2 gc$/(Ki + 2K 2 ), where Cq represents Cq| 5 =o- Therefore the configuration which has a 
larger | eg | will show greater change in energy for a given g. Since the hole configuration in 
Fig. [l|(c) has tetragonal symmetry, which is compatible with the Jahn- Teller distortion, it 
has the largest |cg|. Therefore, as \g\ increases, Fig. |l](c) is more favored than Figs. |l|(a) and 
|l](b). Because the energy changes linearly with g, the phase boundaries are straight lines for 
small g, and have cusps at g = , as shown in Fig. f|. Typical variations of e aa corresponding 
to changing \g\ from to 0.4 are about 2 %. The results indicate that the strain generated 
by substrates can change ordered hole configuration and ordering temperature. 

Our results indicate that the interaction of the electronic state and the lattice can be 
the origin of the charge ordering in this material, even though the details of the results are 
dependent on specific choice of K\, K 2 , A, and A. 

A similar calculation has been done for i?ei/2^4fci/2Mn03, hole concentration 1/2. We 
choose the three hole configurations in Fig. |5| to compare the minimum energies. Each 
configuration has an alternating hole distribution in different set of directions : x, y, and z 
directions for Fig. |^(a), x and y directions for Fig. |](b), and y direction for Fig. ||(c). Figure 
|5](b) is the experimentally observed hole configuration]! As we have done for x = 1/8, we 
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consider both the case where hole and orbital state have the same unit cell, and the case 
where the orbital state is composed of the two hole sublattices. Calculations for A/K\ = 
0.0002, X/Ki = 0.045, and K 2 /K 1 = 0.5 show that when (3 is large, configuration in Fig. 
||(a) is the most favored and Fig. ||(c) is the least favored. As [3 is decreased, the favored 
configuration is changed between (3=0.5 and 0.7. After that Fig. |5|(c) is the most favored and 
Fig. ||(a) is the least favored. When (3 is large, the holes prefer to distribute evenly because 
of the same reason as in x = 1/8 case. In contrast, when (3 is small, electron sites prefer to 
have more neighboring electron sites to gain orbital energy. Our results are not consistent 
with the experimental results for Lai^Ca^MnOa,! which indicate that configuration Fig. 
|5](b) is the ground state. This inconsistency may arise because our model involves only 
localized electrons, while for x = 1/2 the charge ordering state arises from a metallic phase. 
Modifications of our model to include hole hopping are desirable. 

In summary, we have shown that the lattice effect could play an important role in the 
charge ordering transition observed in perovskite manganites. 
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APPENDIX: 



To find the minimum energy we transform E tot in Eq. (|12"D into k space, using Eqs. ([13]) 
(|TED. This leads to the following energy expressions in k space: 



E^/{NK X ) = £ 5\M^ k + 44u 5 + u\L^ k + u£u s + u\P^ k + e\Plu 



k U k 

k 



-^£(1-^)008 6^, (Al) 

where 

4 = K>4^)> ( A3 ) 
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1+ffl 
+f(l- 



3^(1 - cost 



cos k x ) 
><j. cos ky 



-Ea.(i . 



■V) 

cos k x cos k z ) 



sin k y sin fc^. 



|^ sin fc 2 sin /c^ 



sin fc x sin fc^ 



cos k„ 



+J^(1 — COS fcy COS fc 2 ) 

+ — cos /c y cos k x ) 



sin /c 2 sin k y 



\ 



sin fc x sin k z 



sin fcy sin k z 



l + |i(l_ cos ^) 

+J^(l — COS k z COS fc x ) 
+ -j^(l — COS /c 2 COS j 



(-\(l + e ik *) 

-|(l + e^) |, (A6) 

-i(l + e ifcz 

1 - e lk y j . (A7) 









V 





A 


(I 


2K X 


\ 




and N is the total number of Mn sites. We obtain Eqs. flTq) — fl22|) by minimizing the above 
expression with respect to all and w|. Without the second neighbor elastic energy term, 
8? and ii£ minimizing Eq. ( |AT| ) become singular when any of k x , k y , and k z is zero. With 

nonzero K3 this singularity has been uniquely solved for k 7^ 0, while at k — 0, it is not. 

To find the energy term with /c = 0, we take the k — > limit. That corresponds to the 
uniform strain energy, i.e., the energy related to the change of the lattice parameters from 
the original cubic structure. Here the problem of the choice of the limiting process arises, 
because the calculation shows that the different directions of the limiting process of k — > 
have given different energies and different uniform strains. Because the lower energy state is 
favored after all, the appropriate limiting process will be the one which gives the minimum 
uniform strain energy, and it determines the uniform strain also. When K3/K1 -C 1, this 
appropriate limiting process has been found to satisfy the condition of k x , k y , and k z 7^ 0. 
.A.S f £ix cis k% 1 ky ; ciiici k z are nonzero, the limits are different only in the order of K 3 /Ki. 
Therefore, in the K 3 /Ki — > limit, any k — > limit process satisfying the above condition 
gives the correct expression of the minimum energy term with k = 0. It also gives a unique 
uniform strain. 
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FIGURES 



FIG. 1. The three hole ordering patterns for Re 7 / 8 Ak 1 / s MnC>3 considered in our calculations. 
Solid circles represent Mn 3+ , and open circles Mn 4+ . 

FIG. 2. Spring constants: K\ between the nearest neighbor Mn-O, K 2 between the first 
neighbor Mn-Mn, and K3 between the second neighbor Mn-Mn. 

FIG. 3. Phase diagram in A/K\ versus (3 plane for Re 7 / 8 Aki/ 7 MnOs. \/K\ = 0.045, 
K 2 /K 1 = 0.5, and g = 0. 

FIG. 4. Phase diagram in g versus (3 plane for Re 7 / s Aki/ 7 MnOs. \jK\ = 0.045, K 2 /K 1 = 0.5, 
and A/Ki = 0.0002. 

FIG. 5. The three hole ordering patterns for Re 1 / 2 Ak 1 / 2 MnOs considered in our calculations. 
Solid circles represent Mn 3+ , and open circles Mn 4+ . 



11 



TABLES 



TABLE I. Coordinates of site i, orbital states, ionic displacements, and uniform strains for 
the minimum energy configuration of Fig. [j](b), when A/K\ = 0.0002, \/K\ = 0.045, K2/K1 = 0.5 
and (3 = 2.5. 



i 


(nf , nH ,nf) 


9j (radian) 


8f - 65 

* fc=0 


sy - 51 

1 k=0 


Sf - (55 

* k=0 


uf — 

* fc=0 


V V 

u- — u- 

1 fc=0 


uf — u z - 


1 


(0,0,0) 


nnlp mfp 











-0.135 


-0.134 


-0.159 


2 


(I 0) 


1.11 





0.007 





0.135 


0.007 


-0.039 


3 




1.97 











0.004 


0.134 


-0.030 


4 


(I I 0) 


0.03 











-0.004 


-0.007 


0.047 


5 


CO I) 


0.09 





-0.005 


-0.049 





-0.023 


0.011 


6 


(1 1) 


0.09 
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Figure 2, K. H. Ahn et al. 




Figure 3, K. H. Ahn et al. 
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Figure 4, K. H. Ahn et al. 
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Figure 5, K. H. Ahn et al. 



